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ABSTRACT 

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explo¬ 
sion models. We investigate constraints on explosive conditions derivable from measured amounts 
of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized ther¬ 
modynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 
58 Ni/( 54 Fe + 56 Ni), which tends to grow with higher neutron excess and with higher entropy. For 
SN 2012ec, a supernova that produced a Ni/Fe ratio of 3.4 ± 1.2 times solar, we find that burning 
of a fuel with neutron excess i) k 6 x 1CU 3 is required. Unless the progenitor metallicity is over 
5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Su¬ 
pernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, 
at least partially, ejected in the explosion. We find that common spherically symmetric models of 
Mzams 13 M 0 stars exploding with a delay time of less than one second (M cut <1.5 M 0 ) are able 
to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such 
as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower 
limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar 
derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and 
neutrino-neutronization obtained in electron-capture models remain the only viable explanation. 
Keywords: supernovae: general, individual (SN 2012ec, Crab, SN 1987A), nuclear reactions, nucle¬ 
osynthesis, abundances 


1. INTRODUCTION 

Core-collapse supernovae (CCSNe) are the explosive 
deaths of massive stars (Mzams ^ 8 M 0 ). As the core 
of the star collapses to a neutron star or a black hole, 
a shock wave ejects the mantle and envelope at high 
velocities. Nucleosynthesis products from hydrostatic 
and explosive burning are dispersed into the interstel¬ 
lar medium, providing a major production channel for 
the metals in the Universe. 

Comparison of spectral models of standard stellar 
evolution and explosion models with observed nebu¬ 
lar CCSN spectra show encouraging agreement (e.g., 
Dessart & Hillier 2011; Jerkstrand et al. 2012, 2014). 
However, the nature of the explosion mechanism re¬ 
mains unclear. One promising scenario is that of delayed 
neutrino-driven explosions (see e.g., Nordhaus et al. 
2010; Janka et al. 2012; Couch 2013; Bruenn et al. 2013, 
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and references therein). An important test for such mod¬ 
els is detailed comparison of explosive nucleosynthesis 
yields with those inferred from observed nebular-phase 
spectra. A new method to determine the ratio of the 
yields of nickel and iron, after radioactive decays, was 
presented by Jerkstrand et al. (2015, J15 hereafter) and 
applied to several CCSNe. Together with literature val¬ 
ues, a picture emerged that several CCSNe show Ni/Fe 
ratios that are far above the solar ratio. It is of inter¬ 
est to consider how such Ni/Fe ratios arise, and whether 
they offer constraints on explosion models. 

The shock front that travels through the star after the 
core has collapsed compresses and heats the overlying 
layers. The innermost layers experience explosive silicon 
burning which produces iron-group elements such as 
56 Ni, which powers much of the electromagnetic display 
of the SN through its decay to 56 Co and 56 Fe. Silicon 
burning also produces other isotopes such as 57 Ni, 
58 Ni, and 44 Ti, in amounts that depend on the detailed 
properties of the progenitor structure and the thermo¬ 
dynamic conditions. Diagnosis of the yields of these 
isotopes can thus provide constraints on the core-collapse 
process (e.g., Hashimoto et al. 1989; Thielemann et al. 
1990; Woosley & Hoffman 1991; Thielemann et al. 
1996; Seitenzahl et al. 2014; Grefenstette et al. 2014; 
Perego et al. 2015). 

Pioneering calculations of explosive silicon burn¬ 
ing were undertaken by Truran et al. (1966) and 
Truran et al. (1967), and the process was expounded 
upon by Woosley et al. (1973). A key property of the 
burning in typical core-collapse environments is that it 
occurs on sufficiently short timescale (< 1 s) that weak 
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reactions have little time to have a significant influ¬ 
ence (Fowler & Hoyle 1964). The proton and neutron 
numbers are therefore preserved, and as most progen¬ 
itor layers have about equal amounts of protons and 
neutrons, the ash is 56 Ni rather than 56 Fe. The neu¬ 
tron content of the fuel can be characterized by the 
neutron excess 77 = (N n — N p ) / (N n + N p ), where N n 
and N p are the number of neutrons and protons, or 
equivalently by the electron to baryon fraction Y e = 
(1 — 77 ) /2. For many regimes the burning passes through 
a phase of quasi-equilibrium which has weak sensitivity 
to the initial composition but a large sensitivity to the 
value of 77 (Bodansky et al. 1968; Woosley & Hoffman 
1992; Hix & Thielemann 1996, 1999; Meyer et al. 1998; 
The et al. 1998; Magkotsios et al. 2010, 2011). There is 
therefore a direct link between observed iron-group yields 
and three fundamental properties of the explosion; the 
temperature, the density, and the neutron excess of the 
fuel. As common ID progenitor models show the neutron 
excess to vary significantly with mass coordinate (e.g., 
Thielemann et al. 1990, 1996; Woosley & Heger 2007, 
WH07 henceforth), this offers a potential method for 
constraining which layers are ejected and which are not. 
Determination of this “mass cut” in turn reveals infor¬ 
mation on the nature of the compact remnant and the 
explosion mechanism. 

In this work we explore what constraints can be derived 
from measured yields of iron and nickel, as are now avail¬ 
able for several CCSNe (see J15 and references therein). 
Under most burning conditions, the iron comes predom¬ 
inantly from 56 Ni, and the nickel comes predominantly 
from 58 Ni. We focus in particular on SN 2012ec, a Type 
IIP SN with a progenitor detection (Maund et al. 2013), 
a well-sampled light curve (Barbarino et al. 2014), and 
detection of stable nickel lines in the nebular spectrum 
(J15). In §2 we investigate constraints on the nucle¬ 
osynthesis obtained from parameterized thermodynamic 
trajectories over an extensive peak temperature - peak 
density plane. In §3 we consider what progenitor layers 
undergo silicon burning in spherically symmetric stellar 
evolution and explosion models. In §4 we discuss the ef¬ 
fects of asymmetries and neutrino-processing, as well as 
implications of the Ni/Fe ratios measured in SN 1987A 
and the Crab, and in §5 we summarize our findings. 

2. PARAMETERIZED NUCLEOSYNTHESIS MODELS 

We explore models for nickel and iron production us¬ 
ing standard, parameterized thermodynamic trajecto¬ 
ries. This single-zone approach assumes that a pass¬ 
ing shock wave heats material to a peak temperature 
T p and compresses it to a peak density p p . The material 
then expands and cools on a prescribed trajectory until 
the temperature and density are reduced to the extent 
that nuclear reactions cease (freeze-out). Here, we use 
adiabatic (constant T 3 /p) expansion trajectories (Arnett 
1971; Woosley et al. 1973) 

— = = -E. m 

dt 3 r dt t 


T(t) = T p exp(—t/3r) p(t) = p p exp(-f/r) (2) 


with a static free-fall timescale for the expanding ejecta 
(Fowler & Hoyle 1964) 

r = (24 TTGpp)- 1 ' 2 » 446 /p^ 2 s (3) 

Magkotsios et al. (2010) demonstrates by comparison 
with trajectories from several core-collapse simulations 
that nucleosynthesis yields are generally accurate to 
within a factor 2 using this treatment. One may also 
use p(t) in the expansion timescale instead of the peak 
density p pi see Magkotsios et al. (2010) for examples of 
such models. 

Using the code described in Magkotsios et al. (2010) 
(see also Magkotsios et al. 2011), we calculate the mass 
fractions of nuclear isotopes produced by nuclear burning 
for different values of T p , p p , and initial electron frac¬ 
tion Y e . We chose peak temperatures and peak densi¬ 
ties spanning the range of 4 x 10 9 ^ T p ^ 10 x 10 9 K, 
10 4 ^ pp ^ 10 1Q g cm -3 . This parameter space covers 
the conditions encountered in most CCSN models that 
produce any significant amounts of Fe or Ni isotopes. 
The parameter space is sampled with 121 logaritmically 
spaced points, so for any value of Y e we compute the final 
nucleosynthesis at 121x121 points in the (Tp,p p ) plane. 

Our initial composition is the mixture of 28 Si and neu¬ 
trons that give the specified Y e \ this is achieved by mass 
fractions X( 28 Si) = 1 — 77 and X(n) = 77 . The choice 
of initial composition is not important for vast regions 
of the chosen thermodynamic parameter space (using 
e.g. 28 Si and 29 Si to set Y e gives the same results) be¬ 
cause memory of the initial composition is quickly erased 
(Magkotsios et al. 2010). 

For the Y e parameter, we focus attention on Y e = 
0.490, 0.495, 0.497, and 6.499 as these values are repre¬ 
sentative of different shells in the pre-supernova progeni¬ 
tor structure (see Sect. 3). Although there are deep-lying 
layers in the progenitor with Y e less than 0.490, this value 
marks the lower limit below which 56 Ni is no longer the 
major nucleus produced (e.g. Thielemann et al. 1990), 
and is therefore of limited interest for this study where 
we explore SNe with significant 56 Ni production. No can¬ 
didate fuel for silicon burning has Y e > 0.499, and higher 
values than this are not explored. 

2 .1. Overview Of Results 

Figure 1 shows the resulting Ni/Fe mass ratio (after 
radioactive decays) in the ( T pi p p ) plane for Y e = 0.490, 
0.495, 0.497, and 0.499. Contour lines are drawn at 
Ni/Fe = 0.13 and 0.27, bracketing the range derived for 
SN 2012ec in J15 (compare with the solar ratio of 0.056, 
Lodders 2003). For Ni/Fe ratios matching SN2012ec, 
the nickel yield is dominated by 58 Ni, and the iron yield 
is dominated by 56 Ni + 54 Fe, with 56 Ni usually the more 
abundant. We can therefore reach an understanding of 
the observed Ni/Fe ratio by looking at the behavior of 
these three isotopes. 

Figure 2 shows the mass fractions of 56 Ni, 54 Fe, 58 Ni, 
and 4 He, for the case Y e = 0.490 (other Y e values give 
qualitatively similar trends). There are three main pro¬ 
duction regimes. At Tg < 4.5 is the incomplete burn¬ 
ing regime where little iron-group production occurs. At 
Tg > 4.5, log Pp > 7 is the complete burning regime with 
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Figure 1. The Ni/Fe mass ratio (after radioactive decays) in the ( T p ,p p ) plane, for an exponential expansion starting with y e =0.490 
(top left), Y e = 0.495 (top right), 17=0.497 (bottom left), and Y7=0.499 (bottom right). On the x-axis Tg = T p / 10 9 K. Black contour 
lines delineate the value range for SN 2012ec (0.13 - 0.27). Also plotted in the upper left panel are locations of different mass elements in 
multidimensional explosion simulations of SN 1987A (yellow, Wongwathanarat et al. 2010) and of a 27 Mg star (pink, Muller et al. 2012). 


normal freeze-out 9 . In this regime burning is well de¬ 
scribed by a single Quasi-Static Equilibrium (QSE) clus¬ 
ter, in which abundances are largely determined by nu¬ 
clear Q-values. At Tg > 4.5, log p p < 7 is the complete 
burning regime with a-rich freeze-out, where high pro¬ 
duction of ct-particles depletes iron-group yields. Burn¬ 
ing is now described by two separate QSE clusters, as 
the triple-alpha reaction has fallen out of equilibrium. 
The transition region between the normal and a-rich 
freeze-out regimes is called the chasm region, here ma¬ 
terial passes from being dominated by one QSE cluster 
to two. At low densities (logp p < 6) and high ternper- 

9 We adopt here a helium mass fraction X t 1 lie) Rf 0.1 as the di¬ 
viding line between “normal” and “alpha-rich” freeze-out regimes. 


atures (Tg > 7), Y e increases during the burning due 
to weak interactions and mostly protons and a parti¬ 
cles are produced (the so called ap regime). For a more 
refined description of the different burning regimes see 
Magkotsios et al. (2010). 

Figure 3 shows the radiation entropy log (S-y/R), where 
S 7 = AaT 3 / (3 p) and R is the gas constant. The entropy 
values will be of use for later discussion. 

Figure 2 allows a quantitative description of the 56 Ni, 
54 Fe, and 58 Ni yields, and from them the Ni/Fe ~ 
58 Ni/( 54 Fe + 56 Ni) ratio. The neutron excess of 56 Ni is 
zero, whereas for 54 Fe and 58 Ni it is 0.037 and 0.034, 
respectively. Thus, as long as these three isotopes dom¬ 
inate the composition, the neutron excess of the matter 
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Figure 2. The mass fractions of 56 Ni (top left), 54 Fe (top right), 58 Ni (bottom left), and 4 He (bottom right), for an initial Y e = 0.490. 


must be stored in some combination of 54 Fe and 58 Ni. 
Their mass fractions X will then be constrained by 

r) = 0.03 7 X( 54 Fe) + 0.034 AT( 58 Ni) . (4) 

The maximum mass fractions of these two isotopes are 
thus AT( 54 Fe) < 77/0.037 and AT( 58 Ni) < p/0.034. The 
peak plateau values of 54 Fe and 58 Ni production in Fig. 2 
correspond to the mass fractions that give the correct 
neutron excess when that isotope dominates the compo¬ 
sition. 

The neutron excess is predominantly stored in 54 Fe at 
low entropy and in 58 Ni at high entropy. In the nor¬ 
mal freeze-out regime, the Ni/Fe ratio must therefore in¬ 
crease with increasing entropy as 54 Fe is replaced by 58 Ni. 
When the entropy is large enough, both 56 Ni and 58 Ni are 
replaced by a-particles. For high temperatures, weak in¬ 
teractions de-neutronize the matter and a decrease in the 


Ni/Fe yield follows in the ap regime. At lower temper¬ 
ature, where the original neutron excess is maintained, 
the Ni/Fe ratio stays high as 56 Ni is destroyed more effi¬ 
ciently than 58 Ni in strong a-rich freeze-out. 

2 .2. Y e = 0.490 

The Ni/Fe ratio at Y e = 0.490 is shown in Fig. 1 (top 
left). In the complete burning regime, there is a mini¬ 
mum for the ratio at entropy log (S-y/R) « —2; at lower 
entropy the 58 Ni yield increases due to electron captures 
which lower Y e , whereas at higher entropies it increases 
as 58 Ni replaces 54 Fe as the main neutron excess stor¬ 
age nucleus. This growth continues well into the a-rich 
freeze-out regime (log(S' 7 /i?) > +0.5), but is quenched 
in the ap-regime. The maximum 56 Ni production is only 
~40% as 54 Fe and 58 Ni are produced in large amounts 
at these relatively large neutron-rich compositions. 
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Figure 3. Radiation entropy, log S-y/R. 

The band bounded by the contour lines for SN 2012ec 
lies between the normal freeze-out regime and the chasm 
between the normal and a-rich freeze-out regimes (see 
Fig. 4 of Magkotsios et al. 2010). Entrance and exit into 
the allowed band (along increasing entropy) is driven by 
switchover between 54 Fe and 58 Ni, while 56 Ni changes 
little. Using Eq. 4 with r] = 0.02 gives X( 54 Fe) < 0.54 
and X( 58 Ni) < 0.59. The large plateau region of max¬ 
imum 58 Ni production thus gives a Ni/Fe ratio > 1 
(X( 56 Ni) < 0.4), over 18 times the solar value. 

In the incomplete burning regime (Tg < 5), the allowed 
strip between the Ni/Fe = 0.13 and 0.27 contours is very 
narrow and it is unlikely that the burning in SN 2012ec 
occurred precisely under these conditions. The same con¬ 
clusion can be drawn for the narrow strip delineating the 
transition into the ap-regime. 

2.3. Y e = 0.495 

At Y e > 0.495 (r) < 0.01), the neutron excess is too 
small for 54 Fe (p = 0.037) and 58 Ni (17 = 0.034) to dom¬ 
inate the composition, and 56 Ni (77 = 0) is the most 
abundant iron-group nucleus produced. The Y e =0.495 
calculation is shown in Fig. 1 , top right. The Ni/Fe 
band allowed by the SN 2012ec abundance determina¬ 
tions moves to slightly higher entropies relative to the 
Y e =0.490 case. Note how the Ni/Fe = 0.27 contour is 
beginning to move away from the ap-regime. The width 
of this lower density band is not significant at Y e = 0.495, 
but will continue to widen and become significant as Y e 
increases. The plateau region of maximum 58 Ni produc¬ 
tion has X( 58 Ni) = 0.29, giving about six times solar 
Ni/Fe. 


2.4. Y e = 0.497 

The Y e =0.497 calculation is shown in Fig. 1, bottom 
left. This Y e is interesting because it provides a very large 
region in thermodynamic space where the inferred Ni/Fe 
ratio of SN 2012ec is produced. The allowed Ni/Fe band 
covers entropies — 1 < log ( S 1 /R) < + 1 , and crosses into 
the a-rich freeze-out regime. A transition into the al¬ 
lowed band from low to high entropy occurs again as 
the 58 Ni abundance increases when the entropy increases 
past log (S-y/R) ~ —1 (normal freeze-out). However, 
the 58 Ni abundance now levels off to its maximum value 


X( 58 Ni) = 0.18 (from Eq. 4) before the upper bound¬ 
ary (Ni/Fe = 0.27) of the regime is crossed, and gives 
a Ni/Fe ratio that lies within the tolerance interval for 
SN 2012ec, explaning the large size of the allowed region 
at this Y e . 

A qualitative difference to lower Y e values is that now 
a part of the allowed band lies in the a-rich freeze-out 
regime. In particular, the upper Ni/Fe boundary (0.27) is 
crossed (with increasing entropy) not as 58 Ni increases by 
too much (as at lower Y e ), but by a more rapid depletion 
of 56 Ni compared to 58 Ni in the a-rich freeze-out regime. 

A 58 Ni/ 56 Ni ratio of 0.2 ratio corresponds to an elec¬ 
tron fraction Y e =0.497 if the abundances of other nuclei 
are neglegible. If the freeze-out composition is domi¬ 
nated by these two isotopes, then y e =0.497 must be the 
electron fraction of the fuel (assuming Y e stays constant 
during the burning). If other isotopes are present at 
freeze-out, but 56 Ni still dominates the mass fraction, 
Y e = 0.497 represents the maximum allowed electron 
fraction, as no significant amounts of proton-rich iso¬ 
topes (Ye > 0.5) are produced and the other isotopes 
must contribute zero or positive neutron excess. How¬ 
ever, incomplete burning or a strong a-rich freeze-out 
may allow also larger Y e values as 56 Ni is then not the 
main nucleus produced. 


2.5. Y e = 0.499 

The Y e =0.499 calculation is shown in Fig. 1, bottom 
right. There is no regime at normal freeze-out that pro¬ 
duces enough 58 Ni to reproduce the SN 2012ec Ni/Fe ra¬ 
tio, as the maximum 58 Ni fraction X( 58 Ni) = 0.06 (from 
Eq. 4), and thus Ni/Fe < 0.06 (as long as 56 Ni dom¬ 
inates). The only allowed band is restricted to a nar¬ 
row region in the a-rich freeze-out regime with densi¬ 
ties logpp < 5.5. The allowed band has an entropy of 
log(5 7 /i?) = 1.5 — 1.7. At these peak initial conditions 
a-particles dominate the final composition, and X( 56 Ni) 
< 0.3. The allowed band is reached as 56 Ni is depleted 
more strongly than 58 Ni in the a-rich freeze-out. 


2 . 6 . Relation to progenitor density and shock velocity 

Our Y e sensitivity study suggests two fundamentally 
different ways that a Ni/Fe ratio of ^0.2 can be achieved. 
The first is a normal freeze-out burning at low entropies 
(—1 < log (Sy/R) < 0) of a high neutron excess fuel 
(Y e = 0.490 — 0.497). The second is an a-rich freeze- 
out burning at high entropies (log(5 7 /i?) > 0 ) of a 
lower neutron excess fuel ( Y e = 0.497 — 0.499). At nor¬ 
mal freeze-out, the Ni/Fe ratio grows with entropy be¬ 
cause the 58 Ni abundance grows (at the expense of 54 Fe), 
whereas in an a-rich freeze-out the Ni/Fe ratio grows 
because 56 Ni is depleted more efficiently than 58 Ni with 
increasing entropy. 

In the non-relativistic, strong shock limit for a 
radiation-dominated gas (adiabatic index 7 = 4/3), the 
post-shock conditions are related to pre-shock conditions 
by (e.g., Sedov 1959; Chevalier 1976) 


P post — 


Ppre — ^ Ppre 


, 7 - 1 , 

Tpost = 4300 plitv^ 2 K , 


( 5 ) 
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Table 1 

Constraints on p pre and v s in SN 2012ec from application of the 
strong shock conditions. 


Ye 

log ppre (g cm -3 ) 

log u s (km s 1 ) 

0.490 

7.2 - 9.2 

3.2 - 3.6 

0.495 

7.1 - 8.7 

3.4 - 3.6 

0.497 

5.2 - 8.2 

3.6 - 4.6 

0.499 

3.2 - 5.0 

5.0 - 5.5 


where v s is the shock speed. Thus, a [T pos t, Ppost] pair 
maps onto a [/5p re ,u s ] pair, as 

_ 1 

Ppre — y P post 

0 ' w ) cms ' 1 - (6) 

By equating the post-shock conditions with the peak 
conditions in the nucleosynthesis models (i.e., T post = T p 
and Ppost = Pp)-, a determined Ni/Fe ratio constrains 
the allowed ranges of p pre and v s ; the values for SN 
2012ec are listed in Table 1. More neutron-rich (lower 
Y e ) material requires higher pre-SN densities and slower 
shock speeds to achieve the same Ni/Fe production. Note 
that the Ye=0.499 case requires shock speeds exceed¬ 
ing c/3. Electron capture SNe are capable of achiev¬ 
ing shock speeds around this mark as the shock ac¬ 
celerates down their extremely steep density gradients 
(Janka et al. 2008). But SN 2012ec, or any other CCSN 
producing M( 56 Ni) 0.01 M@, must arise from a more 

massive progenitor, and the cores of these are shallower 
(approximately p oc r -3 ), which allows no significant 
shock acceleration (Matzner & McKee 1999). The initial 

velocity scale of the shock is limited to v s ~ yj2E/M < 
0.1 c from the release of 8.8 MeV per baryon in fusion 
reactions (see Scheck et al. 2006, for a discussion of the 
energy budget of the explosion). In a self-regulating ex¬ 
plosion mechanism, such as the neutrino-driven one, the 
shock can also never significantly exceed the escape ve¬ 
locity as the energy deposition shuts off as soon as that 
happens. This gives a constraint v s < \J2GMJr < 0.2c 
for M = 1.4 M 01 and R = 150 km (typical shock radius 
prior to explosion). We conclude that the Y e = 0.499 case 
has no realistic shock scenario associated with it, and can 
be ruled out as a viable scenario for silicon burning in SN 
2012 ec. 

The Y e < 0.497 scenarios all require physically reason¬ 
able progenitor densities and shock speeds. We consider, 
however, that Y e < 0.495 is unlikely because the required 
entropies are over an order of magnitude lower than en¬ 
countered in typical core-collapse explosion simulations 
(a few examples of temperatures and densities obtained 
in core-collapse simulations are plotted in the upper left 
panel of Fig. 1). The much larger regime in thermody¬ 
namic space at Y e = 0.497 instead overlaps with typical 
explosion entropies. We therefore conclude that the most 
natural scenario in which the Ni/Fe ratio of SN 2012ec is 
explained is that the SN burnt and ejected a progenitor 
layer with Y e ss 0.497. 

3. STELLAR EVOLUTION AND EXPLOSION MODELS 

So far our analysis has been focused on single-zone 
models with parameterized thermodynamic trajectories, 


allowing us to explore how the Ni/Fe production ratio 
depends on burning conditions. In this section, we ex¬ 
amine stellar evolution and explosion models from the lit¬ 
erature to explore which progenitors and explosions can 
give the required thermodynamic conditions and neutron 
excess of the fuel derived in Sect. 2. This is a challeng¬ 
ing step since our understanding of both progenitor evo¬ 
lution and explosion mechanisms is far from complete, 
but illustrates the role of the Ni/Fe ratio as an impor¬ 
tant constraint for the testing of both current and future 
supernova models. We limit ourselves in this section to 
spherically symmetric models, providing some discussion 
of multidimensional effects in Sect. 4. We begin by de¬ 
scribing the evolution of the neutron excess during the 
pre-SN evolution. 

3.1. Neutron excess during pre-SN evolution 

Figure 4 shows the final Y e profile prior to core¬ 
collapse of a solar metallicity, 201 isotope, -Mzams=15 
Mq model calculated with MESA, public release ver¬ 
sion 7315 (Paxton et al. 2011, 2013). The bound¬ 
aries of the He, O, Si, and Fe shells are la¬ 
beled. Core helium burning increases rj from its ini¬ 
tial value close to zero ( Y e = 0.5) through the se¬ 
quence 14 N(a, 7 ) 18 F(/ 3 +) 18 0 (a, 7 )(e.g., Couch & Arnett 
1972; Arnett & Thielemann 1985; Thielemann & Arnett 
1985). The 18 0 is burnt further, but in a 77 -preserving 
way. The final neutron excess of the helium burning 
ashes is therefore the neutron excess of 18 O (rj = 0 . 11 ) 
times the mass fraction of material burnt to 18 O (and 
further), which equals 18/14 times the mass fraction of 
CNO, as the CNO cycle converts most CNO to 14 N, and 
the burning of 14 N to 18 0 occurs with close to 100% ef¬ 
ficiency in the core. The CNO mass fraction is about 
2/3 of the total metal mass fraction at solar metallic¬ 
ity (Asplund et al. 2009). Assuming this to hold also for 
other metallicities, we obtain 

^He-burn ash ~ 1® x 014 ) ' ^ 

Core carbon burning further increases 77 
through 12 C(p, 7 ) 13 N(/ 3 + ) 13 C(a, n) 16 0 as well as 
20 Ne(?i, 7 ) 21 Ne(p, 7 ) 22 Na(/ 3 + ) 22 Ne(a, n) 25 Mg(p, 7 ) 
26 Al(/l + ) 26 Mg (Thielemann & Arnett 1985). The in¬ 
crease in 77 is larger for lower-mass cores; about a factor 
of two for a solar metallicity 4 Mq He core but less 
for more massive ones. Note that carbon burning will 
produce a non-zero neutron excess even if Z = 0. At zero 
metallicity, Woosley et al. (2002) find 77 = 1.2 x 10 -3 
and 6.8 x 10 -4 for 15 M 0 and 25 M 0 progenitors, 
respectively. As an approximate formula we may take 

^C-burn ash ~ 10 q 9 ^ 4 ^ ' (®) 

The neutron excess shows no significant changes dur¬ 
ing core or shell neon burning. Core oxygen burning 
increases the neutron excess to 77 = 0.01 — 0.03 (Y e = 
0.485 — 0.495) due to electron capture reactions on 
33 S, 35 C1, 31 P, and 32 S, with higher values for lower- 
mass cores (Thielemann & Arnett 1985). However, 
layers later exposed to explosive silicon burning will, 
for most progenitors and mass cuts, have experienced 
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Figure 4. Y e profile at core-collapse of a solar metallicity 15 M@ 
MESA model. The Fe, Si, and O composition layers are labelled, 
and open circles mark the locations of grid points in the model. 

shell oxygen burning rather than core oxygen burn¬ 
ing. The lower density and higher temperature in 
shell oxygen burning makes electron captures less ef¬ 
ficient and ty is smaller (Thielemann & Arnett 1985; 
Thielemann et al. 1990; Woosley et al. 2002) (see also 
Fig. 4). Thielemann et al. (1990) obtain ry 0.01 
(Y e ss 0.495) in this layer, whereas WH07 have rj ~ 
4 x 10 -3 (Y e « 0.498). In out MESA simultion we obtain 
fj = 6x 10 _3 (Y e = 0.497) in this layer. During core sili¬ 
con burning a large number of electron capture reactions 
increase ry to values ry > 0.06, (Y e < 0.47, see Fig. 4). 
As for oxygen burning, shell silicon burning also gives 
lower neutron excess values, with Thielemann & Arnett 
(1985) obtaining ry = 0.04 (Y e = 0.48) in a 8 M 0 He core, 
similar to the value obtained in the outer Fe core in our 
MESA calculation for a 15 M 0 star (Y e = 0.486, Fig. 4). 

Three distinct zones with very different ?y are thus 
present in the progenitor; the iron core (rj > 0.04), 
the silicon shell (rj ~ 0.01), and the oxygen shell (rj ~ 
10~ 3 x (1 + 21/0.014)). The inferred Ni/Fe ratio in SN 
2012ec is produced for ?y «6x 10 -3 (Sect. 2). There are 
then two candidate origin sites; 1) The silicon-shell 2) 
The oxygen-shell in a very high-metallicity progenitor, 
Z > 0.07. The second scenario requires five times the 
solar metallicity, much higher than the measured metal¬ 
licity of Z = 0.014 — 0.025 (1-1.8 times Zq) in the SN 
2012ec region (Ramya et al. 2007). We conclude that 
SN 2012ec, and likely most other SNe 10 that produce a 
Ni/Fe ratio of several times solar, burnt and ejected part 
of the silicon layer. 

3.2. Explosion 

Explosions in spherical symmetry can be modelled 
using a thermal bomb (e.g., Thielemann et al. 1996), 
a piston-driven event (e.g., Woosley & Weaver 1995), 
or a kinetic bomb (e.g., ChiefH & Limongi 2013). 
Young & Fryer (2007) discuss the differences in the yields 
between thermal bomb and piston-driven methods. In 
the specific case of a neutrino-driven mechanism, the ex- 

10 Even the highest abundances in the inner few kiloparsecs of 
spiral galaxies such as M101 are at most twice solar (Bresolin 2007) 
and a value of 5 times solar would be unprecedented. 


plosions can be modelled with self-consistently calculated 
neutrino luminosities (e.g., Kitaura et al. 2006) or with 
tuneable ones (e.g. Perego et al. 2015). 

The explosion mechanism will take some time fdelay 
after initiation of collapse to send the shock wave off. If 
we specifically define this time as the time at which the 
shock wave has expanded enough to give a post-shock 
temperature below the silicon-burning limit, and denote 
the radius of the shock at that point as -Rsi-burn, the 
total mass that will have been burnt is 

—burn 

M Si —burn = / p{r, £ d elay)47rr 2 dr. (9) 

Jo 

A quite accurate estimate of the size i?si-burn of the 
region that experiences complete silicon burning can be 
obtained from the equation (Woosley 1988) 

= ■^^Si-burn a ^Si-burn > (10) 

Solving for i?Si-burn with Tsi-burn = 5 X 10 9 K gives 

{ E \ 

-Rsi-burn = 3700 km (^ 10 51 erg J • ( n ) 

For prompt explosions the relevant density profile is that 
of the progenitor before any infall, p(r, fdeiay = 0)- If 
the explosion is delayed, then matter has time to accrete 
and the amount of mass inside -Rsi-burn will be higher. 
For example, for Mzams = 13,15,20, and 25 M 0 pro¬ 
genitors, Thielemann et al. (1996, henceforth T96) find 
A?Si-burn = 1-42,1.46,1.70, and 1.79 M 0 for E= 1B (1 B 
= 1 Bethe = 10 51 erg) thermal bomb explosion and a 
delay time of zero. For a delay time of one second, 
these masses increase to -Msi-bum = 1.50,1.53,1.77, and 
1.91 M 0 . 

The mass cut M cut (the dividing point between mat¬ 
ter that falls onto the compact remnant and matter 
that is ejected) cannot be computed ab-initio from 
thermal bomb or piston-driven models. Nevertheless, 
the ejection of even small amounts of iron core ma¬ 
terial results in significantly non-solar abundance pat¬ 
terns (Arnett 1996), and one can therefore argue that 
most SNe should have their mass cuts above the iron 
core (which equals the mass coordinate of the inner 
edge of the silicon shell). This mass coordinate is 
A?Fe—core = 1-18,1.28,1.40,1.61 M 0 for Mzams = 
13,15,20, 25 M 0 in the T96 models. 

For a given explosion energy E and density profile at 
the time of shock passage p(r, tdeiay) we can thus approx¬ 
imate Msi_burn using Eqs. 9 and 11. For a given ejected 
56 Ni mass M 56 Ni,ejected the mass cut M cut is then 

-M C ut = -Mgj_burn dl56Ni,ejected i (12) 

where we have assumed that the mass of the silicon burn¬ 
ing ashes is dominated by 56 Ni. Table 2 lists the resulting 
mass cuts for different progenitors from the T96 model 
grid (using M 56 Ni, ejected = 0.03 M 0 as determined for 
SN 2012ec), and the Y e value in the region between M cut 
and Msi-bum (which becomes the ejected material that 
has experienced silicon-burning). 

Figure 5 (top) shows the relevant quantities from the 
T96 grid, also plotting M cut for non-zero delay times. For 










Table 2 

From the T96 models, the outer mass coordinate for the silicon 
shell Msi.ont (taken as the point where Y e crosses 0.4985), the 
inferred mass cuts M cu t (assuming zero delay-time, E= IB, and a 
Si-burning ash mass of 0.03 Mg), and the Y e values in the 
silicon-burning region for those mass cuts. 


Mzams 

Msi.out 

M cu t 

Y e 

13 

1.49 

1.39 

0.491 

15 

1.37 

1.43 

0.499 

20 

1.63 

1.67 

0.499 

25 

1.64 

1.79 

0.499 



Mzams (M 0 ) 



Mzams (M 0 ) 

Figure 5. The mass cuts for E= 1 B and Af( 56 Ni)=0.03 Mg, com¬ 
pared to the outer edge of the silicon shell Mgi,out (black, solid), 
from the T96 models (top) and the WH07 models (bottom). In 
the T96 models, ejection of silicon-shell fuel occurs only for the 
Mzams = 13 Mq progenitor. In the WH07 models, silicon-shell 
material is not ejected for any progenitor mass. 

E > 1 B, only the combination of a Mzams = 13 M 0 pro¬ 
genitor and a delay time of less than one second ejects 
any silicon-layer material. The silicon-shell Y e in this 
particular model (Y e = 0.491) is too low for SN 2012ec, 
but the value of Y e in the silicon shell depends on de¬ 
tails in the stellar evolution model. For larger Mzams 
and/or longer delay times, the mass cut lies above the 
outer edge of the silicon shell, so only oxygen-shell ma¬ 
terial with Y e = 0.499 is ejected. 

Figure 5 (bottom) shows the same quantities for the 
WH07 models. These models have slightly denser cores 



Figure 6 . The position of M cu t for different values of E, using 
the WH07 models. 

for a given Mzams- The trends are, however, similar. 
The WH07 models show the same behavior of the outer 
edge of the silicon shell approaching the mass cut for 
low Mzams ; in this grid they converge around Mzams = 
12M 0 . The mass cut lies just outside the silicon shell 
here rather than just inside as in the T96 models, but 
small changes in delay time and/or explosion energy will 
move the boundary. 

Using either the T96 or WH07 model grids gives a con¬ 
sistent picture that ejecting the silicon layer at Mzams > 
13 M 0 (and synthesizing M( 56 Ni) ~ 6.03 Mq) would 
require an explosion energy smaller than 1 B. Figure 6 
shows how Msi_burn varies if E is smaller than 1 B. At 
E = 0.5 B, ejection of silicon-layers could occur up to 
Mzams = 15 M 0 , but at larger Mzams still only oxygen 
layers are ejected. Going to E = 0.25 B does allow for 
ejection of silicon-layers at higher Mzams- However, the 
expansion velocities scale as V oc y/E /M e j ecta , and for 
such a low explosion energy we would expect at least a 
factor two more narrow lines than usual. In SN 2012ec 
both photospheric and nebular lines appear as broad as 
in other Type IIP SNe (Barbarino et al. 2014, J15). This 
excludes the scenario of a low-energy explosion of a high- 
mass progenitor, which otherwise may be able to eject 
silicon-layer material (note that we have not adressed 
constraints from the dynamics here). An upper limit to 
the explosion energy for SN 2012ec can also be deduced; if 
E was significantly greater than 1 B no progenitor could 
eject the silicon layer and still produce as little 56 Ni as 
0.03 Mq. 

3.3. Model grids 

The next step is to see what Ni/Fe ratios are pro¬ 
duced in full explosion simulations. We let M(Ni) and 
M(Fe) denote the mass of all Ni and Fe, respectively, 
one year after explosion. SN 2012ec has M(Ni)/M(Fe) = 
0.13 - 0.27 and M( 56 Ni) = 0.02 - 0.04 M 0 (J15). Fig¬ 
ure 7 shows the measured position of SN 2012ec in the 
M( 56 Ni)-M(Ni)/M(Fe) plane compared to the explo¬ 
sion simulations by Woosley & Weaver (1995), T96, and 
Limongi & Chieffi (2003). 

The mass of ejected nickel in these models depends 
on the choice of mass cut/piston mass coordinate, 
which is set manually. Most of the models in the 
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Figure 7. Position of SN 2012ec (M( 56 Ni) = 0.02 - 0.04 M 0 , 
M(Ni)/M(Fe) = 0.13 - 0.27) in the M( 56 Ni) - M(Ni)/M(Fe) 
plane (shaded region), compared to the model grids presented 
by Woosley & Weaver (1995) (WW95), Thielemann et al. (1996) 
(T96), and Limongi & Chieffi (2003) (L03). Also shown is the lo¬ 
cation of SN 1987A. 
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Woosley & Weaver (1995) and T96 grids eject more 
56 Ni than produced in SN 2012ec (so the mass cut is set 
“too deep”). They also produce a lower Ni/Fe ratio. 
Some of the models in the Limongi & Chieffi (2003) grid 
produce lower amounts of 56 Ni, but none has the right 
combination of 56 Ni mass and Ni/Fe ratio. 

The T96 grid uses deeper mass cuts than is appropriate 
for SN 2012ec. For example, the mass cut used for their 
20 M 0 model is 1.61 M 0 , which adds 0.02 M 0 of silicon- 
shell material with r/ = 0.494 to the ejected ashes. This 
gives a dramatically higher Ni/Fe ratio than if the mass 
cut would have been placed at 1.67 M 0 , as needed for an 
ejected 56 Ni mass of 0.03 M 0 . 

For the 13 M 0 model the Ni/Fe value is 0.93. This is 
burning of silicon-shell material with Y e = 0.491. Com¬ 
paring with Fig. 1, the entropy must have been higher 
than what would have been needed to make a ratio con¬ 
sistent with SN 2012ec. Indeed, inspection of our MESA 
simulation shows that the progenitor density at the outer 
edge of the silicon shell is p pie = few x 10 6 g cm -3 so 
the post-shock density (which is about 7 times larger) 
is of order p post = 10 7 g cm -3 . For T ~ 5 x 10 9 K 
the entropies will be higher than needed for Y e = 0.490 
in Fig. 1. For the 25 M 0 star, the Ni/Fe ratio is 0.04, 
too low for SN 2012ec. The burning occurred at lower 
entropy than needed for Y e = 0.499. 

The WW95 explosions use a piston rather than a ther¬ 
mal bomb. With this method, the choice of mass coordi¬ 
nate for the piston combined with a choice of explosion 
energy determines the mass cut (which equals the piston 
mass coordinate plus fallback). In the WW95 models the 
piston is placed at the edge of the iron core. 

Figure 8 shows the Ni/Fe ratio in units of the solar 
value versus M z ams in the WW95 grid. The mass cut is 
inside the silicon-shell for Mzams < 15 M 0 , raising the 
Ni/Fe production to several times the solar value. This is 
not offset by subsolar production at higher Mzams; the 
burning occurs on the large plateau region in Fig. 1. The 
Mzams = 11 — 13 M 0 models give high enough Ni/Fe to 
match SN 2012ec, but the total mass of the ejected ashes 
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Figure 8. The Ni/Fe ratio (relative to the solar value) in the 
WW95 model grid, for E = 1.2 B and solar metallicity. 


is somewhat too large (Fig. 7). The best fitting model 
is 12A, which makes 0.04 M 0 of 56 Ni and gives a Ni/Fe 
ratio of 0.18, in close agreement with SN 2012ec. 

The L03 models are piston-driven explosions and nu¬ 
cleosynthesis yields for different piston motions at each 
Mzams are presented. Most models produce a Ni/Fe ra¬ 
tio around solar. The only models with Ni/Fe approach¬ 
ing SN 2012ec (30B and 35D) only eject trace amounts 
of iron-group nuclei. 


4. DISCUSSION 
4.1. Multidimensional effects 

The ordered onion-like structure of iron-silicon-oxygen 
shells obtained in one-dimensional stellar evolution mod¬ 
els may not be a good approximation of the true 
structure. Multidimensional simulations of the ad¬ 
vanced stages of burning suggest strong convective 
overturns and a rearrangement of the various burn¬ 
ing ashes (Bazan & Arnett 1998; Meakin & Arnett 2006; 
Arnett & Meakin 2011; Meakin et al. 2011; Viallet et al. 
2013; Couch et al. 2015). This opens up a possibility 
that neutron-rich material gets mixed out and resides at 
a larger mass coordinate than in ID models, providing an 
alternative means of making a large Ni/Fe ratio. Inves¬ 
tigation of the explosive silicon burning process in such 
multidimensional progenitors would be of significant in¬ 
terest. 

Also the explosion may involve asymmetries. Direct 
evidence for asymmetries in core-collapse explosions 
are available from imaging of Cas A (Fesen et al. 
2006; Hwang & Laming 2012; Isensee et al. 2010; 
Grefenstette et al. 2014) and SN 1987A (Wang et al. 
2002; Kjaer et al. 2010; Larsson et al. 2013; Boggs et al. 
2015). If the explosion is asymmetric, the entropy for a 
fixed explosion energy becomes higher in the direction 
of the explosion, and the outer mass coordinate for 
complete silicon burning increases in this direction. At 
the same time, we expect deeper layers to be more easily 
ejected in the direction where the energy is focused. 
Thus, one expects the expelled silicon burning ashes 
to originate from a more extended mass range in the 
direction of the explosion. 

Nagataki et al. (1997) studied the influence of 
asymmetries on the 58 Ni/ 56 Ni production using 
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two-dimensional hydrodynamic simulations (see also 
Nagataki et al. 1998; Nagataki 2000). The study 
employed two different mass cuts, the first spherically 
symmetric, and the second sorting test particles by 
energy and making the mass cut at the energy giving 
the specified amount of 56 Ni. For the second (more 
realistic) case, higher degree of explosion asymmetry 
gave deeper mass cuts in the explosion direction, as 
expected. Silicon-layer material can thus more easily be 
ejected in the explosion direction, thereby achieving a 
high Ni/Fe ratio. 

Of significance for the analysis here is also their Case 
B model in which Y e was artificially set to Y e = 0.499 
everywhere. This model illustrates the sole effect of a 
higher entropy achieved in asymmetric explosions (5/ 
increases by a factor of 2-3 from the spherical symmetric 
case to their most asymmetric model). From Fig. 1 we 
expect a quite modest change in Ni/Fe for an entropy 
change of that order; indeed the ratio varies by less than 
30% in the Nagataki et al. (1997) simulations. It is of 
interest to note that the ratio they obtained decreases 
with higher degree of asymmetry. Considering Fig. 1 
again, this is likely because the burning region spans both 
the a-rich freeze-out regime (where Ni/Fe increases) and 
the ap-regime (where Ni/Fe decreases), and the net effect 
is a small decrease. This, in fact, adds further weight 
to the argument that the burning in SNe giving several 
times solar Ni/Fe ratios cannot have occured at as high 
Y e as 0.499. 

4.2. Neutrino-processed ejecta 

Potentially, the material with high neutron excess 
could also come from the hot, neutrino-heated bub¬ 
ble deeper in the supernova core (whose neutron-to- 
proton ratio has been reset by n(^ e , e~)p, p(y e , e + )n, and 
p(e~,u e )n reactions) instead of being directly ejected af¬ 
ter undergoing explosive burning. However, the similar¬ 
ity of v e and v e luminosities and mean energies seen in 
modern simulations generally drives Y e in the innermost 
ejecta above 0.5 (Frohlich et al. 2006). Neutron-rich con¬ 
ditions can only be maintained in some of the early 
neutrino-heated ejecta if they expand rapidly enough, as 
in the electron capture supernova model of Wanajo et al. 
( 2011 ). 

For SNe producing M( 56 Ni) » 0.01 M 0 , this process 
does not provide a likely explanation, however. The to¬ 
tal mass of iron group material originating from the fast 
ejection of neutron-rich, neutrino-processed material is 
limited to a fraction of the mass of the gain region at 
shock revival in this scenario, i.e. to a few times 1O~ 3 M 0 
(as in Wanajo et al. 2011). Moreover, the fast ejection of 
the early neutrino-heated ejecta in electron capture su¬ 
pernovae depends crucially on the special density struc¬ 
ture of their super-AGB progenitors, and it is doubtful 
whether sufficiently short expansion time-scales could be 
reached in more massive progenitors. 

4.3. Nickel isotopes in the Crab 

The Crab has a measured Ni/Fe ratio of 60 — 75 times 
solar (MacAlpine et al. 1989). This extreme ratio can 
be produced in the electron capture scenario described 
above, and indeed nucleosynthesis and expansion dy¬ 
namics are consistent with such an origin (Nomoto et al. 


1982; Kitaura et al. 2006; Wanajo et al. 2009). It is of 
interest to ask whether such a process is a unique solu¬ 
tion, or whether explosive burning without any strong 
neutrino processing of the ejecta may also explain this 
value. 

Inspection of Fig. 1 shows that for Y e > 0.495 (the 
lowest value encountered outside the iron core), the only 
burning conditions that give such a Ni/Fe ratio he in the 
high entropy regime that is not matched by any physi¬ 
cal shock speeds (Sect. 2.6). The most neutron-rich fuel 
Y e = 0.490 can give such a Ni/Fe ratio at lower and more 
typical explosion entropies (log(5 7 /i?) ~ +0.5), but as 
discussed in Sect. 3.1, such a Y e is too low for the silicon 
shell, and would be part of the iron core. It is difficult to 
achieve mass cuts that deep, and many other isotopes are 
produced in extremely non-solar proportions. It there¬ 
fore appears a contrived scenario to explain the Ni/Fe 
ratio in the Crab without the neutrino processing occur¬ 
ring in an electron capture event. 

4.4. Nickel isotopes in SN1987A 

In J15 it was demonstrated that taking the optical 
depth of [Ni II] 6.636 pm into account, the Ni/Fe ratio in 
SN 1987A derived from nebular-phase line luminosities is 
around solar, suggesting no or small ejection of silicon- 
shell material. This is consistent with the analysis in 
Kozma & Fransson (1998), where a model with a nickel 
mass of 1.3 times solar relative to iron gave reasonable 
fits for [Ni II] 6.636 pm and [Ni II] 10.68 pm. These re¬ 
sults supersede initial estimates of Ni/Fe <0.5 times so¬ 
lar based on analytical formulae assuming optically thin 
emission (Rank et al. 1988; Wooden et al. 1993). The 
abundance of 57 Ni could be directly determined from 
gamma-ray lines, giving a 57 Ni/ 56 Ni ratio of 1—2 times 
the solar 57 Fe/ 56 Fe ratio 11 (Kurfess et al. 1992). A sim¬ 
ilar range was inferred from infrared lines (Varani et al. 
1990) and light curve models (Fransson & Kozma 1993, 
2002; Seitenzahl et al. 2014). 

Figure 9 shows the production of Ni/Fe (in units of 
solar) and 57 Ni/ 56 Ni (in units of solar 57 Fe/ 56 Fe) for the 
Y e = 0.499 case from the parameterized thermodynamic 
trajectories of §2. A Ni/Fe ratio of 0.5 — 1.5 times solar 
and a 57 Ni/ 56 Ni ratio of 1—2 times solar is reproduced for 
burning at post-shock densities of p p = 10 6 —10 8 g cm -3 , 
which are typical values encountered in simulations (Fig. 
1). The revision of the Ni/Fe ratio in J15 due to the 
calculated optical depths in the [Ni II] 6.636 pm line is 
therefore important in allowing a consistent solution for 
the production of nickel isotopes in SN 1987A. 

A lower Y e moves the fit region to higher peak densi¬ 
ties. For example, Y e = 0.497 would require p p > 10 8 
g cm -3 to produce the SN 1987A values (not plotted, 
but see also Woosley & Hoffman (1991)). This density is 
larger than obtained in any published explosion mod¬ 
els. Thus, the most consistent picture is that explo¬ 
sive burning in SN 1987A occurred in the oxygen shell, 
and that the mass cut therefore was outside the sili¬ 
con shell. For a Mzams = 20 M 0 spherically symmet¬ 
ric progenitor model this is at ~1.65 M 0 in the T96 
models and at ~1.76M 0 in the WH07 models. From 
this we expect the minimum value of the mass cut to be 

11 ( 57 Fe/ 56 Fe) 0 =O.O23 (Lodders 2003) 
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Figure 9. The Ni/Fe ratio in units of the solar ratio (top) and 
the 57 Ni/ 56 Ni ratio in units of the solar 57 Fe/ 56 Fe ratio (bottom), 
from the Y e = 0.499 simulation. 

around 1.7 Mq. The baryonic neutron star mass was con¬ 
strained to 1.2—1.7 Mq from the neutrino burst (Burrows 
1988). This is marginally consistent with the constraints 
from the silicon burning, and disfavors a progenitor mass 
larger than 20 Mq (but allows for smaller). 

5. SUMMARY 

We have analyzed the production of nickel and iron 
(after radioactive decays) in explosive silicon burning, 
and the constraints on progenitor structure and burn¬ 
ing conditions derivable from measured yields. We have 
focussed on the implications of Ni/Fe ratios of several 
times the solar value, as recently reported for the Type 
IIP SN2012ec (J15), and previously for the Type Ic SN 
SN 2006aj (Maeda et al. 2007; Mazzali et al. 2007). 

Using a nucleosynthesis network with parameterized 
thermodynamic trajectories, we have computed the 
nickel and iron yields (after decays) as function of peak 
temperature, peak density, and neutron excess rj of the 
fuel. The Ni/Fe ratio is for most regimes dominated 
by the production ratio of 58 Ni/( 54 Fe + 56 Ni). In the 
normal freeze-out regime, this ratio grows with increas¬ 


ing entropy as the neutron excess storage switches from 
54 Fe to 58 Ni. When this switch is complete, the ra¬ 
tio reaches a plateau value of Ni/Fe = 77/0.034; fuels 
with 77 = 0.034 x Ni/Fe provide a large region of al¬ 
lowed thermodynamic conditions. Smaller values of 77 
can also achieve the production, but only in a limited re¬ 
gion of strong a-rich freeze-out where the 58 Ni/ 56 Ni ra¬ 
tio grows due to efficient depletion of 56 Ni. We find that 
this process is not likely to be responsible for produc¬ 
ing Ni/Fe ratios several times solar as the necessary en¬ 
tropies would require unphysical shock speeds. Higher 
values of 77 can achieve the production, but only in re¬ 
stricted regions of low entropy that are not encountered 
in typical explosion simulations. We conclude that to 
produce a Ni/Fe ratio of a few times solar, burning at 
normal freeze-out, or moderate a-rich freeze-out (helium 
mass fraction X( 4 He) < 0.2), of a fuel with 77 ~ 6 x 10 -3 
(Y e ~ 0.497) is required. 

The derived neutron excess value can be linked to 
the location in the star where explosive silicon burn¬ 
ing occurred. Models for the progenitor structure show 
that the neutron excess successively increases inwards 
from the oxygen shell (?) ~ 2 x 10 “ 3 ) to the iron core 
(i)>4x 10~ 2 ). Only in the silicon shell is the neutron 
excess in the range needed for a Ni/Fe production of a 
few times the solar value (77 ~ 6 x 1CU 3 ). An excep¬ 
tion is if the metallicity of the star exceeds five times 
solar, in which case the neutron excess in the oxygen 
shell increases to 77 ~ 6 x 10~ 3 . Such a metallicity 
can be ruled out for SN 2012ec from HII region spec¬ 
troscopy (Ramya et al. 2007). SN 2012ec, and other CC- 
SNe producing large Ni/Fe ratios, therefore likely burnt 
and ejected part of their silicon layers, which is a key 
constraint for explosion models. 

In spherical symmetry, a given progenitor structure, 
explosion energy, delay time, and measured amount of 
56 Ni, defines the position of the mass cut. Lower-mass 
stars have relatively thick silicon shells that more easily 
encompass the mass cut. We find that Mzams ^ 13 
Mq progenitors exploding with a delay time of less than 
one second are able to eject part of their silicon layers. 
Such a progenitor for SN 2012ec is in agreement with the 
estimate of Mzams from modelling of the hydrostatic 
burning ashes (Mzams = 13 — 15 M 0 , J15). Higher- 
mass progenitors only eject oxygen shell material with 
moderate a-rich freeze-out, giving a Ni/Fe ratio close to 
the solar value. In particular, the measured amounts of 
56 Ni, 57 Ni, and 58 Ni in SN 1987A appear consistent with 
the burning and ejection of pure oxygen-shell material. 
This in turn translates to a lower limit for the mass cut 
of this SN of 1.7 Mq. 

Asymmetry in the explosion can qualitatively change 
the Ni/Fe ratio by two means; by more easily ejecting 
deeper-lying silicon layers in the direction of the ex¬ 
plosion, and by achieving a stronger a-rich freeze-out. 
Published simulations (Nagataki et al. 1997) show that 
the entropy effect is insufficient, but the first mechanism 
can achieve a high Ni/Fe ratio. An asymmetric explo¬ 
sion is a plausible explanation for the high Ni/Fe yield 
in SN2006aj, an X-ray flash SN (Mazzali et al. 2006; 
Maeda et al. 2007; Mazzali et al. 2007). 

For the Crab, a very high Ni/Fe ratio of 60-75 yimes 
solar has been reported (MacAlpine et al. 1989). We 
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find that this extreme value cannot be reproduced 
by any realistic-entropy burning outside the iron core, 
and neutrino-neutronization obtained in electron-capture 
models (Wanajo et al. 2009) remains the only viable ex¬ 
planation. 

In conclusion, it is clear that constraints on both 
progenitor structures and explosion dynamics can be 
obtained from silicon-burning yields determined from 
nebular-phase spectra. The exact location of the sili¬ 
con burning (i.e. the mass cut) has strong ramifications 
for iron-group yields obtained in galactic chemical evolu¬ 
tion models, and the constraints derived from SN 2012ec 
and several other CCSNe provide important input for 
such modeling. Further observations and modelling of 
nebular-phase SNe, combined with nucleosynthesis mod¬ 
elling, will enable us to make further progress in under¬ 
standing the elusive SN explosion mechanism and the 
origin of the iron-group elements. 
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